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I. INTRODUCTION 



Moyal star-product formalism, more generally deformation quantization, provides useful 
tools for describing the phase space formulation of quantum mechanics jl|, Q]. Spin and 
relativistic quantum mechanics can also be considered in this framework by making use of 
an additional fermionic star product [3] . It is seen that a combination of the Clifford product 
as fermionic part and Moyal star-product as bosonic part plays an important role in these 
regards. This approach can also be adopted to realize the methods of supersymmetric 
,uantu,n mechanic Qlflg on a pha.e =pace. Hen^elde. ha. already .hown how .uch^ 
approach can be accomplished on a two-dimensional {2D) phase space in a recent study [8|. 
The main goal of the present work is to show that such a programme can be realized on a 
4D phase space in its full generality and rich applicability potential. 

In this paper, we first study the deformation of the complex Clifford algebra defined 
on a 4D phase space of a classical Hamiltonian system. The deformation is performed in 
composing the components of Clifford forms with the Moyal star-product. This enables us 
to obtain two iso-spectral matrix Hamiltonians, each being a sum of factorized terms. The 
Hamiltonians depend on four real- valued phase-space functions and they are general enough 
for modeling various systems. Moreover, the Hamiltonians are also doubly intertwined 
in such a way that the intertwining matrix-valued phase-space functions are divisors of 
zero in the resulting deformed algebra. By composing the intertwining matrices with their 



conjugates, we obtain a constanto. 
type models of quantum optics 



motion for each system. Well-known Jaynes-Cummings- 



10, 



which are the main subject of recent experiments 



12], emerge as special cases 



to understand the quantum nature of atom-field interactions 
of our method. These are studied in detail as the illustrative examples. Their phase-space 
characteristics such as their spectra, eigenspinors and related Wigner functions as well as 
their classical limits and constants of motion are obtained in an autonomous way, without 
any reference to the standard tools of quantum mechanics. 

The paper is structured into two main parts. The first part consists of the following 
three sections where our method is fully developed. In Sec. II the Clifford algebra structure 
of AD classical phase space and its complexification are studied and the basics of Moyal 
star-product are briefly reviewed. The deformation of the algebra and its application for 
establishment of supersymmetry techniques on the phase space are taken up in Sec. Ill 
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and its physical consequences are studied in Sec. IV. Sec. V, which constitutes the second 
part of the study, is entirely devoted to the applications of the method in quantum optics. 
Concluding remarks are given in the final section. 

II. CLIFFORD ALGEBRA ON A PHASE-SPACE 

Any vector space on which a non-degenerate inner product is defined has a Clifford 
algebra structure. Tangent (or cotangent) space of any (pseudo) Riemannian manifold is 
the best place where efficiency of Clifford calculus can be observed [l3^. In this regard the 
symplectic manifold that we shall consider in this and the following section can be any 
even- dimensional manifold. However, in order to be able to do quantization consistently in 
such a framework, one needs the star product of deformation quantization (see subsection 
IIC) which requires globally defined coordinates. For this reason and for physically relevant 
applications, we shall confine the investigation to AD flat cases. Our method can easily be 
extended to an arbitrary flat symplectic manifold. For definiteness, we take our phase space 
to be a AD symplectic manifold A4 topologically equivalent to and denote the linear 
space of complex-valued, smooth (differentiable to all orders) functions defined on it by JF. 
Real- valued elements of JF are the classical observables of Hamiltonian mechanics. 

At each point of A4 two AD vector spaces, the tangent space and (its dual) the cotangent 
space are naturally defined. Elements of the latter are 1-forms acting linearly on the ele- 
ments of the former. Ai carries two distinguished structures defined by two non-degenerate, 
second rank tensor fields: the symplectic 2-form Q (which is also closed) and the Euclidean 
symmetric metric tensor g. They are non-degenerate in the sense that each defines a vector 
space isomorphism between the tangent and cotangent spaces. The elements associated by 
these isomorphisms are called symplectic and metric dual of each other. The metric dual of 
a vector field x is a 1-form x defined, for a given vector field y, as x{y) = iy{x) = g{x,y). 
Here iy denotes the so-called interior derivation (or interior product) whose action on an 
arbitrary /c-form (a totally anti-symmetric covariant tensor) (3 is defined, for arbitrary vector 
fields Xi, . . . Xk-i, by {iyl3){xi, . . . , Xk~i) = kj3{y, xi, . . . Xk-i)- ^ and g at each point of a 
given neighborhood endow the tangent space with structures of a symplectic vector space 
and an inner product space, respectively. According to the Darboux theorem jl^, in each 
neighborhood of Ai one can define the canonical coordinates (q, p) = {(li,(l2,Pi,P2) which 
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lead to f2 = dqjAdpj. Here A denotes the bilinear and associative exterior product which 
satisfies 

a A (3= {-ly'f^Aa, 
for arbitrary j-form a and /c-form (3. Q eventually leads to the Poisson bracket (PB) 

2 

[F, G]p = Y,i^.Fd,.G - d,,Fd,fi) (1) 

k=l 

of Hamiltonian mechanics. We henceforth use the abbreviation dq^ = d/dqi and denote the 
generic elements of T by capital Latin letters. On the right-hand side of (1), the following 
ordinary commutative and associative pointwise product of functions is essential: 

(FiF2)(q,p) = (F2Fi)(q,p) =Fi(q,p)F2(q,p) =F2(q,p)Fi(q,p) . 



A. Clifford Bundle and Complexification 

The linear space of all forms constitutes the exterior algebra with respect to exterior 
product, which (like the tensor algebra) is a Z-graded associative algebra. Although, unlike 
the tensor algebra, the exterior algebra is finite dimensional (2"D when the dimension of 
M. is n) its Z-gradation is inherited from the tensor algebra. In particular, any fc-form for 
k > n is zero and the zero element of the exterior algebra is the only homogeneous element 

n 

of every degree greater than n (see also [13[ pp. 7-8). 

The Clifford product *ci defined for a 1-form x and an arbitrary form j3 by 

x*cl3 = xAP + i:,l3 , (2) 

turns the exterior algebra into a Z2-graded associative algebra called the Clifford algebra 



(also known as the geometric algebra 
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16| . See also 17| for applications in Hamiltonian 



mechanics). By associativity, the rule (2) suffices to completely determine *c on arbitrary 
forms. Z2-gradation means that the whole algebra IS, clS Qb linear space, a direct sum of 
the spaces of the odd and even forms such that the latter has a subalgebra structure. The 
Clifford product of a j-form a and fc-form (3 is, in general, an inhomogeneous form consisting 
of a sum of £-forms such that i = j + j + k — 2^ . . . — k\. The exterior bundle (union of 
exterior algebras) equipped with the *c'-product in the fibres is called the Clifford bundle. 
The Clifford commutator, defined by 

[a,(3]c = a*c P - P *c a , 



is bilinear, antisymmetric and satisfies the Jacobi identity. Each fibre of the Chfford bundle 
acquires a Lie algebra structure with this bracket. In the context of Clifford algebra, we shall 
mainly adopt the conventions of [3], except for the representation of the Clifford product. 
More conventional versions of this product are x *c (3 = x V (3 and juxtaposing the factors. 

Tangent spaces can be complexified by replacing the field of real numbers M by the field 
of complex numbers C. Then g can be extended by C-linearity to the complex-valued, 
symmetric and non-degenerate (guaranteed by the non-degeneracy of g) bilinear map g'^. 
As C is algebraically closed, g"^ is not characterized by any signature (even if g had any 
non-Euclidean signature) and hence the structure of complex algebra depends only on the 
dimension. In our case, we consider 2*^D real Clifford algebra C4^o(IR) which is isomorphic 



to the algebra of 2 x 2 quaternion matrices 



13| . p. 80). But from here on, we shall deal with 



its complexification C4(C) which is known to be isomorphic to the algebra of 4 x 4 complex 
matrices. 

B. The Primary SUSY Structure 

The 2^D complex algebra C4(C) is generated by 1 and by the complex orthonormal basis 
e^, e^, e^} such that 

*c + *c = '^5^'" , (3) 

where the Kronecker symbols denotes the components of the (inverse) metric g'^{e\ e^) = 5^^. 
We start out our analysis by the Clifford 1-form 

uj = W^e^ + W2e^ + Pie^ + Pse^ , (4) 

whose components are real-valued phase-space functions 

W, = iy,(q, p), P, = P,(q, p) , J = 1, 2. 

A rather general Hamiltonian function in two dimensions can now be written, in terms of 
uj, as a Clifford product 

H=\^*c^ = \{Pl + Pi) + \{W,^ + W2^) . (5) 
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We next introduce the 1-forms 

g ^ ^{e^ + te'^) , g = ^{e^-ie^). 
These are all nilpotent of order 2 



and satisfy 



^ 9*c9 , 


f *C f ^ ^g*cg : 


(6) 


{/, f}c = 


2 = {9,9}c , 




{9j}c = 




(7) 


{9,f}c = 


= {9: f}c , 





where {, }c denotes the Clifford anti-commutator. The relations (6) and (7) imply that the 
set {2~^/^/, 2~^/^g', 2~^/^/, 2~^/^^} constitutes a Witt basis of the complexified cotangent 
space. The first two and the last two elements of this basis span two isotropic subspaces 
(where g'^ induces the zero bilinear map) whose direct sum is the whole cotangent space. 
Then in terms of the complex- valued functions 

C,^^{W^ + lP,), C2 = ^iW2 + iP2), (8) 

and their complex conjugates Ci,C2, we define 

q^^Cif + C2g, q+^Cj + C2g, (9) 

such that (4) can be rewritten as u; = g+ -|- g_. By using (6) and (7), one can easily verify 
that q± are also nilpotent, and together with H, they close in a simple supersymmetric 
algebra structure 

q±*cq± = , 

H = ^{q_, q+}c , (10) 
[q±,H]c = 0. 

The Chfford algebra structure of the exterior bundle of A4 enabled us to see the super- 
symmetry (SUSY) structure in the corresponding classical system. Such systems possessing 
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the fermionic (anti-commuting elements) degrees of freedom, in addition to the usual bosonic 
ones, are known as pseudoclassical models (see for an extensive list of references). They 
serve as the classical limits of quantum systems having both kinds of degrees of freedom. 
However, the SUSY structure emerging above does not seem to be promising much at this 
stage since the third relation of (10) is trivially satisfied. Indeed, if is a 0-form and therefore 
Clifford commutes with all forms. In general, the third relation of (10) must be a conse- 
quence of the nilpotent nature of so called supercharges q±. In what follows, the algebra 
will be deformed in such a way that the last three relations with a new Hamiltonian, which 
turns out to be an inhomogeneous even Clifford form comprising a 0-form and a 2-form, will 
manifest a genuine SUSY structure. For this purpose, the basics of Moyal star-product and 
its corresponding bracket are briefly reviewed in the following subsection. 

C. Star Product and Moyal Bracket 

In the canonical (q, p) coordinates, the Moyal ^-product on JF is defined by 

1 ^ ^ ^ ^ ^ 
* = exp [2^^^^ ( - )] ' (11) 

where h is the Planck constant and d-, d are acting, respectively, on the left and on the 
right. This product is bilinear, associative, and obeys the relation 



(Fi * F2) = F2 * Fi , (12) 
under complex conjugation. In terms of the *-product the Moyal bracket [, ]m is defined as 

[F,G]m = F*G-G*F , (13) 
for all phase-space functions. Note that in view of (12), we have 



[F,G],j = -[F,G]m. (14) 

In particular, the Moyal bracket of two real-valued functions is a purely imaginary-valued 
function. The most important properties of the ^-product and the MB are the following 
limiting relations 

lim F * G = FG , 

liml[F,G]M = [F,G]p, 
R->o in 



which hold for generic ^-independent phase-space functions. These reveal the fact that the 
associative *-algebra and the Lie algebra structure of JF given by the MB are, respectively, 
deformations (in the sense of Gerstanhaber [isl) of associative algebra structure of T with 
respect to pointwise product and of Lie algebra structure determined with respect to PB. 
All quantum effects are encoded in the ^-product with respect to which the real elements of 
T are promoted to the status of quantum observables. 

III. MOYAL-CLIFFORD ALGEBRA 

So far the components of Clifford forms were commuting quantities since they were multi- 
plied by the ordinary pointwise product. However, one can go over to the non-commutative 
or quantum case by demanding that the coefficients are to be multiplied by the Moyal *- 
product. This can be achieved by combining (2) and (11) together. The resulting associative 
product and algebra will be referred to as the Moyal-Clifford (MC) product, denoted by *mc 
and as the MC-algebra, respectively. We will directly apply this product to the above for- 
mulation. In doing so, we would like the first relation of (10) to remain intact with respect 
to this new product as well. 

A. SUSY structure by deformation 

It is easy to show that 

g+*MC'?+ = {Ci^C-^uf *c9 1 

q_*MCq- = [Cl,C2]Mf *c 9 = -[Ci,C2]Mf *c 9 , 

where in the last equality we have used (14). Thus g+ and g_ are nilpotent with respect to 
the MC-product if and only if 

[Ci,C2]m = 0. (15) 
In terms of W/s and -Pj's, this condition amounts to 

m, - [Pi, P2U = ^m, pi]m - ^[w,, (16) 

Since W/s and Pj's are real-valued, the left hand side of (16) is, in view of (14), purely 
imaginary, while the right hand side is real. So the condition (16) is, in fact, equivalent to 
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the following two conditions: 



[W^,W2]m = [Pi,P2]m, (17) 
[W^,P2]m = [W2,Pi]m. (18) 

We now define the SUSY Hamiltonian as 

Hs = q-}MC , (19) 

which implies that 

[q±,Hs]MC = 0. (20) 

In terms of a;i = a; and UJ2 — —i{q+ — q~) which obey {wi,W2}mc — 0, Hs can be factorized 
as 

Hs = -^CUl *MC = ■^<^2 *MC <^2 ■ (21) 

As will be shown in the following subsection, in an appropriate matrix representation of the 
Clifford algebra, cuj are Hermitian supercharges and the symmetry determined by the above 
super-algebra is 2-extended supersymmetry denoted also hy N — 2 SUSY (where N stands 
for the number of Hermitian supercharges) . 
Using (9), (19), conditions (17), (18) and 

[C, Ci]m = t[W,, P,]m , [C2, C2U = t[W2, P2]m , 

[Ci,C2]m = m,W2]M-i[W2,Pl]M , 
[C2,C,]m = -[Wi,W2]M-i[W2,Pl]M , 

Hs can explicitly be evaluated as 

Hs = H, + ^{[W,,P,]Me'' + [W2,P2W 

+ [W2, PiUie'^ + e'') + [W^, W2]M{e'' + e^)} , (22) 

where H^ should be read as i?*l with 1 being the unit element of the Clifford algebra and 

2H^ — {Ci,Ci}m + {C2,C2}m 

= Pi* Pi + P2* P2 + Wi*Wi + W2*W2 . (23) 



Here {, }m stands for the anti-Moyal bracket and we have adopted the abbreviation e^^ — 
*c which for the product of different orthonormal basis elements becomes e^^ = A e^, 
etc. Note that Hg is an inhomogeneous even Chfford form whose bosonic part is a zero 
form and the remaining fermionic part is a 2— form. The existence of all possible 2- form 
basis elements in the fermionic part is a reflection of its generality. 



B. Matrix realization 

In terms of the 2x2 Pauli matrices 

yioy \^°/ yo-i 

and of the 2x2 unit matrix 1, we shall use 









( 


i(T3 \ 














1 




I -i(73 


/ 



(24) 





i(J2 1 


















\ -ia2 


1 









for the complex Clifford basis defined by (3). This is a representation such that all basis 
matrices are Hermitian, while and e*^ are symmetric and and are antisymmetric. If 
the representations (24) are used in (22), we obtain 

/ i^i \ 

Hs= { =i/i7r+ + i/27r-, (25) 

\0 H2J 

where 7r+ = diag{l, 0) and 7r_ = diag{0, 1) denotes the non-primitive projections and 

Hi = H,l + HiF , H2 = Ha + H2F . (26) 

HjF^s represent the following fermionic parts 

HiF = '-B+as , (27) 
H2F = ^S_(73 - i[W2, PiUai - i[Wi, W2\m<J2 , (28) 
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with 

B± = [W^,P^]m±[W2,P2]m . 

In deriving these relations we made use of cxi £72(73 = il and its consequences. Note that Hi 
is diagonal. 

Having obtained matrix realizations of the Hamiltonians we should emphasize that both of 
the Hi and H2 are Hermitian. Firstly, in view of (14) H^ given by (23) is a real- valued phase- 
space function. Secondly, as and the other two coefficient functions of H^f consisting 
Moyal brackets are purely imaginary valued, the presence of the imaginary unit i in (27) 
and (28) ensures the Hermiticity of the both fermionic parts. 

Equations (15) (or equivalently (17) and (18)), (20) and (22) show the main differences 
between the Clifford algebra and the Moyal- Clifford algebra. These restrictions given by 
(15) arise in order to make q± nilpotent in the deformed case and play important roles in 
the rest of this work. The resulting algebra is a genuine SUSY algebra and the resulting 
Hamiltonians have non-classical parts. The common bosonic part if* of Hi and H2 has 
H given by (5) as its classical limit: lim^j^^o-^^* = H. However, their fermionic parts are 
different and have no classical limits despite the fact that the coefficient functions of the Pauli 
matrices have classical limits: all Moyal brackets in (27) and (28) reduce to PB brackets in 
the classical limit after dividing by ih. 



IV. INTERTWINING, ISO-SPECTRAL PROPERTY AND CONSTANTS OF MO- 
TION 



In terms of 



Li = Ci{iai + a2)-iC2{l-a^)=2i\ ° ° |, (29) 

Ci -C2 

C \ 

L2 = Ci(i(7i + (72)+iC2(l + (73) =2i I ' , (30) 

Ci I 

the matrix representations of the supercharges q± are found, from (9) and (24), to be 

1 Li , 

5+ = ^ ]=<!-'■ (31) 



V2 I -Lo 
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The matrices of uji and uj2 are Hermitian. Nilpotency of g+ implies that 

Ll *MC L2 = = L2 *MC Ll , (32) 

i.e. Ll and L2 are divisors of zero with respect to the *MC-product, which denotes, from 
now on, the star product of matrix-valued phase-space functions. We should note that the 
second equality of (32) results from the matrix product of factors iai + 02 and (1 ± as) 
appearing in (29) and (30), while the condition (15) is essential for the first equality of (32). 
On the other hand, (20) implies the following double intertwining relations: 

L2 *MC Hi = H2 *MC L2 , (33) 
Ll *MC H2 = Hi *MC Ll . (34) 

One can easily verify that [Di *mc -^2)^ = L)\ *mc L>\ holds for arbitrary matrix valued 
functions -Dj's. In view of this fact, the Hermitian conjugates of (33) and (34), or equivalently 
the relation [q-^H^j^ic = 0, gives the following additional intertwining relations 

lI *mc H2 = Hi *Mc L\ , (35) 
L\ *mc Hi = H2 *Aic L\ . (36) 

As is evident from (32), L\ and L\ are also divisors of zero. By virtue of (19) and (31), we 
also obtain 

Hi = ^{Li *MC L\ + lI *mc L2) , (37) 
H2 = -{L{*mc Ll + L2 *Mc lI) , (38) 

which show that each partner Hamiltonian can be written as a sum of factorized terms. 

Let us call a nonzero 2x1 matrix- valued function \1/ a star-eigenspinor of a 2 x 2 matrix- 
valued function T corresponding to the star-eigenvalue A if and only if T *mc ^ = -^^^ This 
definition implies that *mcT'^ = A\E'^, and that \1/ is nonzero if and only if *mc ^ 7^ 0. 
Then the standard theorems for Hermitian operators, in particular the theorems concerning 
the reality of spectra and the orthogonality of eigenfunctions corresponding to different 
eigenvalues, are also valid in the present context. We use the term spinor in the usual 
sense of 2 x 1 matrix eigenfunction and remark that this term has a wider meaning in the 
nomenclature of Clifford algebra. Bearing these facts in mind, we now return to the physical 



implications of the intertwining relations |7|. 
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The first important implication is the fact that Hi and H2 are iso-spectral, that is, they 
have almost the same spectra. More concretely, if is an eigen-spinor of Hi with eigenvalue 
A, then L2*mc'^ and L\*mc'^ are also eigen-spinors of H2 with the same eigenvalue provided 
that is not in the MC-kernel of L2, or L\. In view of (34) and (35) the analogous remark 
is valid for the eigen-spinors of H2. For brevity, we refer to the following section for more 
details and return to another implication of the intertwining relations. 

Multiplying one of the relations (33-36) by Lj (or Lj) from the left and comparing the 
resulting expression with a similarly obtained one from the others, lead to the following 
matrix valued functions: 

Ri = Li *MC L\, 82 = L\ *mc L2 , (39) 
i?2 = L2 *MC lI, Si = L\ *mc Li . (40) 

These commute with Hi and H2 in the following way: 

[Ri,Hi]mc = 0=[S2,Hi]mc, (41) 
[R2,H2]mc = = [Si,H2]mc ■ (42) 

However they are not independent since AHi = R1+S2 and 4if2 = -R2 + 'S'i. That is, if there is 
no explicit time dependence, each system has, together with the Hamiltonian, two constants 
of motion. The explicit forms of the constants of motion for Hi are S2 = 2iJi(l + a-^) and 
Ri = 2ifi(l — (T3). Since Hi is diagonal, 5*2 and Ri are its projected forms. 

The method developed so far contains four real-valued phase-space functions. Apart 
from conditions (17) and (18), there is no constraint on these functions. Therefore, the *- 
products of functions appearing in our formulae contain, in general, countably infinite terms 
each characterized by a positive power of h. This generality enables us to study quantum 
properties of many physically relevant systems on a classical phase space. However, as we are 
about to do in the following section, for physically relevant systems some of these functions 
should be restricted to finite polynomial functions of the canonical coordinates. In such a 
case a comparison with the existing literature of supersymmetric quantum mechanics, in 
which majority of application have been carried out in the usual operator formulation of 
quantum mechanics, can be made. In this regard it is interesting to observe that the forms 
of our Hg, Hi and supercharges are similar, up to some permutations of rows and columns 



of matrices, to that found in 



19| 
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V. APPLICATIONS 



In terms of spin raising and lowering matrices a± — {a\ ± 1(72) /2, H2F given by (28) can 
be rewritten as 

H2F = ^5.(73 + V2{Aa^ + Aa.) , (43) 

where 

A = [W2, Ci]m , A = -[W2, Ci]m . (44) 

These are the phase-space analogues of the bosonic raising and lowering operators when 
is a (real) constant. To give illustrative examples we shall (from here on) restrict 
the study to such cases, take h = 1 and denote the second term of (43) as 

HjciA) = V2{Aa+ + Aa^) , (45) 

when .A]m = 1. In such a case, Hjc is the phase space analog of the well-known, fully 
quantum mechanical Jaynes-Cummings ( JC) Hamiltonian of quantum optics which describes 
a bipartite (a two-level atom and a field) system. In the parlance of quantum optics, a± are 
known as the atomic transition operators and (73 as the inversion operator. The first Aa^ 
term at the right-hand side of (45) represents an absorption process in which destruction 
of a field quanta and a transition to a higher atomic level take place. The second term 
Acr_ describes an emission process in which creation of a field quanta is accompanied by a 
transition to a lower atomic level. 

Let us see how the phase-space version of such a model together with its supersymmetric 
partner naturally emerges from the method developed above. 

A. Example 1: Jaynes-Cummings- type systems 

A particular simple choice of Wj and Pi amenable to various physical applications seems 

to be as 

Wi^P2 , Pi = g2 , = qiP2 - q2P\ ■ (46) 

Then, in view of \Wi, W2\m = Wi and \W2-i Pi\m = —iflii the conditions (17) and (18) imply 
that the most general form of P2 is 

P2 = P1P2 + ?i?2 + , (47) 
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such that Ki = Ki{qi,pi). It is now straightforward to verify, with Ci = {p2 + "1(12) / from 
(44), that 

A = — )={qi + ipi) , [A, A]m = 1 . (48) 
v2 

Hence, A and A are bosonic lowering (annihilation) and raising (creation) phase-space func- 
tions, respectively. One can also define another such bosonic pair by 

B = -^{q2 + ip2) , [B,B]m = 1, (49) 

which Moyal commute with the above pair. W2 represents the angular momentum perpen- 
dicular to the gig2-plane and, together with P2, can be rewritten as 

W2 = -i{AB - AB) , P2 = AB + AB + Ki, (50) 

Let us define the phase-space number functions 

Na = A*A = AA-^, Nb = B*B = BB-]-. (51) 

By making use of 

Fi(q) * F2(q) = Fi(q)F2(q) , Gi(p) * ^^(p) = Gi(p)G2(p) , 
and G *jv/ G = for G = a ■ p + -Fi(q) (a is a constant vector), we also find 

H2F = \B_a^ + Hjc{A) , 

5± = i[~l±2{Nrs-N^ + X{)\ , (52) 
= H,, + 2Ns{Nj, + l) + Y^ , 

where Hj^ = {2N_a + l)/2, and 

Xi = -'-[W2,K,]m, (53) 
Fi = ^{AB + AB, K,}m + * K, . (54) 



H2 is of the type of a JC Hamiltonian j^, [lo|, [ill] describing a two-level atom interacting 



with a quantized, two-mode electromagnetic field of which only the ^-mode directly interacts 
with the atom and hence causing transitions between levels. Hj^ in is the energy of the 
^-mode and the remaining part 
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represents the energy of the other mode as well as the interaction between the modes. The 
super-partners are 



= {H,^t + Hj,)l + {H^ - Nis)as , (55) 
H2 = {H,^t + Hj,)l + {Hs - iV^)a3 + Hjc{A) , (56) 

where His = {'^Nq + l)/2, we have taken Ki = and hence Hint = '^Nq{Nj^ + 1). In such 
a case the interactions between the modes are through the number functions. By choosing 
K_\ such that Xi and Yi are nonzero, one can also model certain types of mode interactions 



20l |. Even in the simplest case when Ki is a nonzero constant, Xi is zero but the Yi term, 
which becomes Ki{AB + AB) + {Kf/2), represents a coherent exchange of photons between 
the modes. In any case Hi has no JC-type atom-field interaction term. 

B. Eigenvalues, Eigenspinors and Wigner Functions 

Hi is diagonal and algebraically depends on the complete set {A^^l, N^l, (T3} of three 
mutually commuting (in the MC sense) matrix-valued phase-space functions. Using the 
shortened notation (z| . . . ) = (q, p| . . . ), their common eigenspinors will be denoted by 

(z|j, n_A, nts) = \j) ® (z|n^, ub) , (57) 

Here n_4, = 0, 1, 2, ... , stand for the numbers of mode-quanta and j{= 1, 2) labels the 
bare states of atom such that |j)'s, like the spin-up and spin-down states of a spin-1/2 
system, satisfy crsll) = — |1),(T3|2) = |2) and 

a_\l) = = a+|2) , 

^-|2) = |1) , 
a+|l) = |2) . 



In (57) the real- valued functions (z|n_4,ng) represent the diagonal Wigner functions 2ll. |22 



(z|n^, ns) = ^^^"^ *M *m (z|0, 0) A^^ *m , (58) 

for two-mode field. Here (z|0,0) denotes the vacuum Wigner function defined by 

A*M (z|0,0) = = B*M (z|0,0) . 
16 



The normalized (its integral all over the phase space set to 1) solution of these equations is 
7r~^ exp — (pi + Pi + q1 + Qi)- Explicit functional forms of the higher level Wigner functions 
are proportional to j23 | 

L„^(4^^)L„g(4Si3)(z|0,0) , 
where Lk stands for the Laguerre polynomials. These are products of the harmonic oscillator 



Wigner functions which are well known 24i\ . 
One can also verify 



A*M {2\nj^,nB) = - , 

A*M {z\nj(,ni3) = ^/nXTT {z\nj^ + l,ns) , (59) 
Na*m (z|n^,nB) = n_A (z|n^,nB) • 

Similar relations hold for Ns, B and B. In view of these relations, the eigenvalues of Hi 

Hi *Mc {z\^,nA,nB) = Ai„^„g(z|l,n^,nf5) , 
-f^i *Afc (z|2,n^,rie) = A2„^„s(z|l, n^, rig) , 

are easily found, from (55) and (59), to be 

-^In^ns = ^^6(2^^ + 3) , A2„^„e = (ng + l)(2n^ + 1) . (60) 

For riA = = riis, we have Aioo = and A200 = 1 which represent the bare energy levels of 
the atom. In this case. Hi is simply (1 + 0-3) /2. The level Xin^o = is infinitely degenerate 
whereas all the other levels are finitely degenerate. Provided that > uq, the level A2n^ns 
is the higher level of the whole atom-field system. These degeneracies arise from the fact 
that we have taken the frequencies of the modes to be equal (degenerate modes). Evidently, 
these discussions can be extended to general case in which the frequencies are different. 
For this example, we have C2 = i\/2BA, Ci = —iB and by virtue of 

2cr+o-_ = 1 + as , 2a_a+ = 1 - (T3 , (61) 

the intertwining (matrix-valued) functions, given by (29) and (30) can be rewritten as 

Li = 2ScT_(l + V2A(J+) , lI = 2B{1 + V2Aa_)a+ , (62) 
L2 = 2B{1 - V2Aa+)a. , lI = 2Ba+{l - V2Aa-) . (63) 
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Now it is easy to verify 

L2*MC {'2'\'^-,rij,,nB) = = *MC (z|2,n^,nB), 

which simply follow from the actions of a± on the atomic bare states. However, one needs 
to be more careful for the verifications of 

L| *mc {z\l, nA,nB) = <l>„_^„g(z) , 

L2*MC {z\2,nA,nB) = 2VriBTT *n^ns(z) , 

where 

^n^nsi^) = V2nB + 2 (z|l,n^ + l,ne-l) + (z|2, 71^,710-1) , (64) 
^n^nslz) = (z|l,n^,7T,B + 1) - A/2n^ (z|2,n^ - + 1). (65) 

As a consistency check, one can directly verify that these are the (unnormalized) eigen- 
spinors of H2 

H2*MC ^TiAnBi'^) = ^InAUB^nAnBi'^) ^ (66) 

-^2 *MC *n^nB(z) = '^2n^nB*n^nB(z) , (67) 

with the same eigenvalues as given by (60). These fully establish the iso-spectral property 
of Hi and H2, which was compactly expressed by the SUSY algebra given by (19-21). 

In the nomenclature of quantum optics the *-eigenstates given by (64) and (65) are the 
phase-space analogues of the so-called stationary dressed states or JC-doublet. The states at 
the right-hand sides of (64) and (65) are known as the bare states of the atom-field system. 
The latter are the phase-space version of the product states of the bare atom and field states. 
Finally, we compute, with the help of (62) and (63), the constants of motion 

Ri = 8Nb{H^ + 1)(7_(7+ , (68) 
Si = A{Nb + l)[(7+(7- + 2NBa-a+ + HjciA)] , (69) 
such that [Ri, Hi]mc = = [Si, H2]mc- 

C. Example 2: Non-resonant JC-type interactions 

For Wi — pi, Pi — —qi and W2 as in (46), the conditions (17) and (18) are satisfied, 
provided that 

-P2 = P1P2 + ?1?2 + K2 . 
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Here K2 — ^2 (52,^2) is an arbitrary real function of its arguments. Hence -\/2Ci — P\ — iqi 
and 

[W2,C,]M = -{Q2-tP2)/V2^B, 

where B and B are defined by (49). As a result, by interchanging the lowering functions of 
subsection V.A such that {A,B) <-> {B,A), or by direct computations, we find 

H2F = '-B^a^ + HjciB) , 
Hjc(B) = V2{Ba+ + Ba^) , 

S± = t[l±2{Ni3-N_A + X2)] , 
^ Hb + 2N^{Nb + 1) + ^2 . 

W2 and the first term of P2 are still given by (50), and X2, Y2 are defined as in (53) and (54) 
provided that Ki is replaced by K2. For K2 — 0, the Hamiltonians and the corresponding 
eigenvalues are 

Hi = {Hint + Hb)1- {Hb - Nj,)a^ , 
H2 = {Hint + H^)1- {Hj, - Ns)as + Hjc{B), 
Ain^ns = {riA + l)(2ne + 1), A2n^„B = nA{'^nB + 1). 

These are related to (60) by the interchange (1, n^, ns) ^ (2, ng, n^). 
Now some remarks for the above applications are in order. 

The JC-type Hamiltonians originate from the interaction of the electric dipole moment 
of an atom with a quantized light in dipole approximation and, in general, both Hjc{A) 
and Hjc{A) take part in such an interaction. To see this we should emphasize that dipole 
moment has non-vanishing matrix elements only between states of opposite parity. There- 
fore, the atomic bare states |1) and |2) are assumed of opposite parity and hence the dipole 
moment is proportional to (7+ -|- (T_ = |1)(2| -|- |2)(1|. On the other hand since a quantized 
single mode cavity field is a multiple of ^ -|- ^ then the dipole energy is proportional to 

{a+ + a-){A + A)^ Hjc{A) + Hjc{A) . 

However, the resonant processes described by Hjc{A) are more efficient than the non- 
resonant processes described by Hjc{A)i especially in the case of a quantized single-mode 
light. Neglect of the latter is usually called the rotating wave approximation. On the other 
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hand, in the case of multi-mode interactions, both types of terms have important physical 
implications. Here we will be content with pointing out that the above two examples can be 
combined and then extended in various ways to study the generalized models of quantum 



optics, such as A^-atom JC models 



25| (also known as Dick models 26j) and spin-boson 



systems [27]. In particular, time can be included and the dynamics of transitions between 
atomic levels for both partner Hamiltonians can be studied in the phase space. 



VI. CONCLUDING REMARKS 



The Moyal- Clifford algebra, that is, Clifford algebra endowed with the star product of 
deformation quantization, provides a unified framework to realize the SUSY techniques in 
a classical phase space, to see relations of various models described by matrix Hamiltonians 
and to study their phase-space characteristics in great detail. To emphasize some merits 
of such an approach we recall that in the usual operator formulation of quantum mechan- 
ics, Wigner functions are indirectly obtained from wave functions through some convolution 
integrals which, apart from simple few cases, are not easy to cope with. But, in the deforma- 
tion quantization they are directly obtained from star-eigenvalue equations. Moreover, by 
realizing the SUSY techniques in this framework Wigner functions of some delicate states 
such as the dressed states can be easily obtained by the intertwining functions. 

In the presented examples the phase space is spanned, instead of the usual canonically 
conjugate coordinates of a particle system, by the conjugate amplitudes of mode functions of 
a quantized electromagnetic field that are dynamically equivalent to conjugate coordinates 
of a mechanical oscillator. The form of Hi and H2 and the presented examples may be 
convincing enough that a lot of physically relevant (charged or uncharged) particle systems 
can be identified as spacial cases. As much as exhaustive search in this direction seems to 
be rather involved and will be deferred to a later publication. 
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